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We propose a simple phenomenological model describing vibrational dynamics in a class of ma- 
terials (like clathrates), which are crystalline compounds with closed cavities in their structures. In 
the spirit of a minimalist approach, our model includes only the minimum of ingredients necessary 
to capture the generic features of this class of compounds as brought to light by experiments. We 
consider only two kinds of particles: strongly coupled host atoms, which form the structural cavities, 
and guest atoms confined in the cavities. In spite of the confinement, for relatively small amplitude 
vibrations of the guest atoms, their interaction with the host particles is assumed to be weaker 
than between the host atoms. We calculate self-consistently the vibrational mode (phonon) line 
broadenings for the empty and filled (by the guest atoms) cavities. We show that the soft phonon 
dispersion sheet, which appears for weakly bonded guest atoms, yields to large mode broadening 
via three waves anharmonic mode coupling. In turn the Umklapp part of this broadening reduces 
drastically the phonon thermoconductivity coefficient. We conclude that for a system with a soft 
mode component in the vibrational spectrum phonon broadening is always larger and thermal con- 
ductivity is always appreciably smaller than for a system without this soft mode (provided all other 
characteristics are the same). Although our model is minimal, in the sense of ignoring the complex- 
ity of actual guest-host structures, when properly interpreted, it can yield quite reasonable values 
for a variety of measurable quantities. 

PACS numbers: 05.60.Cd, 44.10.-|-i,44.05.+e 
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I. INTRODUCTION 



Among the various types of so-called guest - host systems (see e.g., the recent review [l| containing numerous 
relevant references) one of the simplest type can be constructed from two kinds of particles, the strongly coupled host 
atoms forming structural cavities, and guest atoms, which can be captured in the cavities. The term guest atoms 
already indicates that they should not play a strong structural role although actual host structures may become 
unstable in the absence of the guest. We may, therefore consider the situation of guest particles that are confined 
within the cavities, but apart from this, are only weakly interacting with the host atoms. The interest to study these 
systems reflects both, their potential practical importance (related to the possibility of selectively tuning mechanical, 
electrical, or thermal properties), and the associated theoretical challenges. However, complexity of composition of 
existing guest-host systems (like clathrate hydrates, or skutterudites) that underlies the rich physical properties, poses 
a formidable problem to theoretical and computational modeling efforts. Although a number of detailed experimental 
investigations and_sophisticated calculations and numerical computations have been published over the last few years 
(see, e.g., papers [2] - l^, which partially will be commented upon in our paper, and much more can be easily added 
to this list), there is still a clear need for a simple (but yet non-trivial) theoretical model with predictions, which can be 
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FIG. 1: Cartoon ID guest-host structure 



confronted with experimental observations. The fact (perhaps a bit taste-dependent) is that more or less realistic and 
complete theoretical or numerical investigations of existing guest - host systems (with about 50 atoms per elementary 
cell) do not yield simple intuitive (yet based on firm physical principles) explanations of the vibrational excitation 
broadening and the considerable reduction of thermal conductivity upon filling of the host structural cavities with 
guest atoms. In spite of the fact that the basic concepts describing vibrational excitations and thermal transport 
were born long ago and can be found in the textbooks [13], many questions of principle remain to be settled 



for guest-host systems (see e.g. the discussions about Einstein rattlers in the very recent papers or about 

electric dipole force contributions in [^). In part this frustrating situation is just due to the lack of a simple and 
tractable analytical model. Our motivation for adding one more paper to the topic is precisely to propose such a 
model. 

The outline of this article is as follows. In Section |TT] we formulate our model of a simplest guest - host system and 
calculate dispersion curves. In Section [llll we derive self-consistent equations for vibrational mode broadening, which 
are solved numerically in Section [TVl Based on this solution we compute also thermoconductivity coefficients. Section 
fVlis devoted to a discussion, and summary of our main results. 



II. PRELIMINARIES. 



We shall recall certain basic facts about the guest-host systems. A generic structural feature of this class of solids 
is that their lattices are formed by the strongly interacting host atoms or molecules. On the contrary the guest 
molecules can be located only in the cages formed within the host structure. The basic interaction between the guest 
molecule and the host lattice is generally weaker than between the host molecules, i.e. the elastic force controlling 
the guest particle's vibrations is somehow small. However, for larger displacements the interaction between the guest 
and host lattice becomes stronger holding the guest particles firmly within the cage. The guest-host interaction can 
in this sense be assimilated to an inverted hard-core potential, soft in the center and hard outside. The weak elastic 
coupling property leads to the presence of low-frequency or soft dispersion sheets in the system, which correspond to 
the vibrations of the guest molecules within the cage. For the sake of simplicity, let us sketch first a one-dimensional 
cartoon (see Fig. 1) corresponding to our model structure (later on in the next section we will do a proper 3D 
calculation). Essential for our model is that the guest- host structure elementary cell contains two kinds of particles: 
structural atoms H (host) and guest atoms G. In equilibrium the H atoms occupy positions xh — 2n + 1, and there 
are only nearest neighbor interactions between the H atoms. Guest atoms G are located between the atoms H at the 
equilibrium positions xq = 2n and they interact only with the host atoms forming their cages. 

The elastic energy for such a system can be written as 

= i ^ (ko(u2„+1 - U2n-l)^ + '«l(M2n - U2ri-l)^ + '«i(m2„ - U2n+l)^) , (1) 
n 

where w„ is a displacement of atom n, kq and ki are the host -host and guest -host force constants. Soft guest modes 
require that ki < kq. 

For the ease of algebra we are using dimensionless units where both masses {H and G atoms) are set equal to 2. 
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FIG. 2: Spectrum of eigenmodes for the elastic constants kq = ki = 2 (ID model). 



The dynamical equations describing the eigenmodes and eigenvectors of this system read as 



n+l 



K0(2U2„+1 - U2n-1 - U2n+3) + Ki(2m2„+1 - U2n - U2n+2) , 
Kl(2u2n - U2n-1 " U2n+l) • 



(2) 



Using the standard Ansatz for the excitations of the form 



this system of equations reduces to 



U2n+i(i) = fli exp{-iujt + iq{2n + 1)), 
'U-2nit) — 0-2 exp(— iwi + i2qn) 



(2ko sin^(g) — cj^ + Ki)ai — ki cos(q)a2 = , 
— Ki cos((7)ai + (kq — t^^)a2 = . 



Nontrivial solutions exists at uj — wi (q) (acoustic branch) or lu — W2 (q) (optical branch) given by the relation 



kq sm'^{q) ^ k\ cos2((j) + sin''((j) . 



(3) 



(4) 



The spectrum of eigenmodes as defined by equation ^ is presented in Figures O and [3] In Fig. [5] the force constants 
are choose as kq = ki = 2. In Fig. [3] kq = 2 as before but the guest-host force constant is appreciably smaller, 
Ki — 0.5. Calculating the eigenvectors of the system ^ it is easy to show, that the region of small wave-vectors 
within the acoustic branch corresponds to in-phase oscillations of both host and guest atoms. For the case presented 
in Fig. [31 the region of large wave- vectors of the lower "soft" branch corresponds to oscillations of dominantly the 
guest atoms. The flat part of this branch resembles that of Einstein-like independent rattling oscillations, which are 
widely discussed in the literature [l| - Despite a certain similarity we would like to stress that, in our model the 
flatness of the branch reflects the soft dynamics of the guest atoms, and not decoupled dynamics of Einstein oscillators 

m- 

As can be easily demonstrated the property controlling the form of the spectrum ^ is the minimal value A^„j„ of 
the gap A^(q) — w^iq) — w\{q) between the two branches. It reads as 



1/2 



(5) 



(cp. with the direct absolute gap defined as [w^"'^)'^ — (w™*")^ — ki). Therefore, as had to be expected the ratio 
ki/kq is (up to a scaling factor) the only relevant parameter in our model. Upon decreasing the ratio ki/kq, one gets 
a more extended "flat" part of the spectrum (what is supposed to be a "rattler" -like decay channel increasing mode 
broadening). Even more, the reader may be surprised to learn that for sufficiently small ki/kq < 1/7 the large q 
(short wavelength) upper branch group velocity, could become larger than the long wavelength acoustic mode (sound) 
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FIG. 3: Spectrum of eigenmodes for the elastic constants kq = 2, hi — 0.5 (ID model). 

speed. In these unusual conditions one would have phonon broadening of this branch via interactions to the lowest 
lying long wavelength acoustic mode, which for normal downward sloping modes are forbidden due to momentum 
and energy conservation rules. However, upon decreasing the ki/kq ratio, we simultaneously increase the gap ([5]) 
and therefore the mode anharmonic coupling (see the next section [TTTl for more details) to provide comparable with a 
larger value of the ki/kq broadening, should be much stronger. 

Thus for strong broadening there is a certain optimum value of the ratio ki/kq , not too small (to have a reasonably 
strong mode coupling) and not too large (to have a reasonably extended flat spectrum). Instead of a full systematic 
study of these issues we perform in this work a few numerical estimates to narrow down the focus of our problem. In 
what follows just this optimal range of the model parameters will be investigated. 

Summarizing up to this point, for our minimal model in the harmonic approximation we have at our disposition 
one physical parameter ki/kq that allows us to tune two phenomena (namely, mode flatness and mode softness, 
however not independently) that contribute to vibrational broadening and thermoconductivity reduction. One might 
generalize the model to include also guest - guest coupling. However, because it is not a major goal of our work to 
achieve agreement between the results obtained with our phenomenological model and experimental measurements, 
we defer the generalization to future publications. Nevertheless since the present understanding of mechanisms leading 
to thermoconductivity reduction in guest - host systems is incomplete, even our minimal model may be an appropriate 
tool for working out typical trends that in principle can be tested in experiments. 



III. VIBRATION EXCITATION SPECTRUM 



Phonons are crucial for the understanding not only of vibrational properties but equally of a material's electric and 
thermal transport. They are also at the heart of a broad range of interesting electro-mechanical effects. In particular 
for the transport properties phonon decay plays an important role. The problem of phonon decay has a rather long 
history, and the phenomenon has been treated in quite some depth in order to rationalize neutron scattering data 
concerning collective excitations. Theoretical concepts on anharmonic vibrational modes in crystalline materials have 
been developed in the seminal papers by Peierls [3], Pomeranchuk fl7l|. Herring fTs*!, and Klemens [l^ (see also 
early works on the decay of an optical phonon into two acoustic phonons via anharmonicities (20| . [2l| . and a more 
recent publication [l^]). The short version of the formal machinery starts as its initial step to uncouple the coupled 
lattice oscillators, and then to express the Hamiltonian for the vibrating lattice as a sum of separate oscillators. Free 
phonons (in our model with the dispersion laws (|4])) result from the quantization of these oscillators. In the harmonic 
(quadratic over lattice displacements) approximation there is no interaction between phonon states, and an arbitrary 
number of phonons may be present in each mode. The anharmonic terms present the processes involving three and 
more phonons describing their decay or coalescence. To find the broadening for the phonon dispersions calculated in 
the previous section |TT] in the harmonic approximation one has to include phonon-phonon interaction self-consistently. 
The diagram technique is an appropriate tool to treat such a problem. We skip almost all general technical details, 
and focus here on ingredients indispensable for our calculations. 

To be more specific we take into account only three wave interactions, because usually these interactions are 
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FIG. 4: One-loop diagrams for phonon mode broadenings 



responsible for the mode broadening [l7| - [l9|. In the symmetric case these three phonon couplings are vanishing, 
we thus have to introduce asymmetry into our simple model in order to achieve non-zero three-wave interaction. 
The asymmetry can be related to a off-center position of the guest atom within the cage, or merely non symmetric 
spring constants between the guest atom and the cavity walls. In this case the anharmonic higher order couplings 
do not introduce new physics, but only renormalize parameter values of the coupling terms. Since these spring 
constants are not well known, we regard their values as phenomenological coefficients. In this sense our model should 
be treated as a working hypothesis. Comparison of the qualitative predictions resulting from this hypothesis with 
experimental or microscopical numerical observations will show whether and when this hypothesis is justified. Note 
also that it is a priori no problem to include also fourth order anharmonic couplings into the same self-consistent 
treatment. This will not affect our qualitative and semi-quantitative conclusions much (albeit at the expense of rapidly 
increasing complexity). We do believe that with the aim of answering first simple questions transparency is worth 
a few simplifications. Guided by this principle we consider in what follows only three wave couplings and only high 
temperatures T > T^, where Td is the Debye temperature. 

In canonical phonon variables the Hamiltonian of the system can be presented as 

= (wiq)aiq)a*iq)+ J g ^ (a(g)a(pK + p) + c.c.)) . (6) 

The vertex V for the interaction of three phonons with quasi wave- vectors q, p and q -I- p and frequencies w, v and 
cij + // is proportional to the product 

V^^Ek±^V. (7) 

\J LLiv(y: + v) 

For simplicity we suppose that the parameter y is a constant. For the sake of the skeptical reader we should note, 
that if the above mentioned asymmetry of the guest - host interactions or/and equilibrium positions is somewhat 
small, the three-wave interaction vertex V has to be small as well. It is not to say that guest-host vibrational mode 
coupling is small. Indeed, as a guest atom approaches the wall of the cavity, repulsive interactions become dominant 
and the host atoms in the vicinity of the wall are pushed back to decrease the unfavorable interactions. Consequently, 
host-guest vibrations are strongly coupled, even though the direct attractive potential between the guest and the host 
lattice is small. To avoid too many unknown parameters we are not introducing explicitely this asymmetry coefficient 
in the ([6]), ([7]). All the more that in this study we are predominantly concerned with only relative values of phonon 
mode broadenings and thermoconductivities for pure host and guest-host systems. 

The peculiarities of the phonon mode broadenings are related to the existence of the soft mode, and these effects 
are determined by inter-band coupling (i.e., the decay of one optical or acoustic phonon into two acoustic phonons, or 
into one acoustic and one optic phonon, and inverse coalescence processes, satisfying the simple selection rule that the 
created phonon must lie in a higher branch than at least one of the destroyed phonons) integrated over the full range of 
wave vectors. In the one-loop approximation, the broadening of the phonon spectra can be determined from diagrams 
of the type schematically shown in Fig. 01 where the lines describe phonons which decay or coalescence. We take 
into account all possible processes involving three phonons (not only inter-band ones, mentioned above). Calculating 
these diagrams we end with the following system of equations to find the broadening of the phonon spectra 

7l(g)= r^'^V^2/ (7l(p)+7l(9+J5))Ww'l(p))-^K(9+P))) 



j-kI2 27r V(wi(g)+t«i(p)-u'i((?+p))^ + (7i(p)+7i(9+P))^ 
(71 (p) + 71(9 - + n{wx(j))) n(wi(g -p)))/2 



(wi((?) - wi(p) - wx{si~v)Y + (71 (p) +7i('7 -P))^ 

(71 (p) +72(g+p))('^(wi(p)) -r^(w2(g+p))) 
(wi(q) +wi(j)) - W2{q+pW + (71 (P) +l2{q+pW 
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(71 (P) + 7i('7 - P))(l + n{wi{p)) + n{wi{q-p)))/2 



{W2{q) - wi{p) - Wi{q~p)Y + (71 (p) +71(9 -P))^ 

(71 (P) + l2{q - p)){l + n{wi{p)) + n{w2{q - p))) 
{W2{q) - wi{p) - W2{q-p)Y + (71 (p) +l2{q-p)Y 



(8) 



Here subscripts 1 and 2 stands for acoustic and optical branches of the spectrum, n(w) is the Bose factor, n{u}) = 
l/(exp{huj/kB/T) — 1). To stay in contact with experiments we consider only the high temperature case when the 
Bose factor turns into the Boltzmann law n(Lu) ~ ksT / [huj) 3> 1. One more comment is in order here. Of course 
physical phonon damping is determined by the functions 71 and 72 calculated at the actual mode dispersion laws 
(|3|) (so-called "on-shell" phonon broadening). However as an intermediate step of calculations, when we wish to 
include self-consistently the effects of phonon damping, " off-shell" functions 7 can enter the above integral equations 
([5]). We check that our self-consistent equations are stable with respect to off-shell virtual phonons, and thus their 
contributions are irrelevant for our work. 

Then, the equations ([5]) can be rewritten as 

^,(,) = ^ ipf^^)( ii{s){ii{p)+ii{s)){Ms)-w,{p)) 

^ Wi{q)J-^/2 ^ V(wi(q) +u)i(p) -u)i(s))2 + (7i(p) +7i(s))2 

, Pi (s) (71 (p) + 71 (s) ) (u^i (s) + wi b) )/2 



{wi{q) -wi{p) -Wi{s)Y + (71 (p) +7i(s))2 

fi{s){li{p)+l2{s)){w2{s)-w,{p)) 
{wi{q) + wi{p) - W2(s))2 + (71 (p) + 72(5))^- 

V ,2r ,f fUsKli{p)+l2{s)Kw2{s)-w,{p)) 



/ ^ ^ M / ' , .2/ J2[S)[ll(P) +12(S))[W2(S) -Wl(p)) 

™ Mq)J-./2 ^^'^^'y{w2{q)+W,{p)-W2{sW + {jl{p)+j2{sW 

I f? (g) (71 (P) + + {p))/2 

{w2{q) - wi{p) - wi(s))2 + (71 (p) -I- 71 (s))^ 

fHs)i7l{p) + I2{s))iw2{s) + Wljp)) 

{W2{q) -wi{p) - W2(s))2 + (71 (p) +72(5))^ 



(9) 



where A = TV^/{2tt), /i,2(p) = \p\/wi^2{p) and s = q + p. To take into account Umklapp processes it is necessary 
to extend the functions w{q) and 7(9) into the region [— 37r/2 , 37r/2] as a periodical repeat of the first Brillouin zone 
[— 7r/2, 7r/2]. Equations (O give the imaginary part of the one-loop correction to the bare Green functions, which 
determines the decrements. For simplicity we neglect the real part of this correction, which only changes the bare 
dispersion lawsjH and treat the system ^ self-consistently in the functions 71 (q) and 72(5). The system can be solved 
numerically, and as a result we can calculate the broadenings of the two branches (|4]) as a function of parameters A 
and the ratio ki/kq. 

These self-consistent equations ^ are derived yet for the one-dimensional model shown in Fig. 1, and the result 
overestimates the contributions from small wave-vectors (although it is worth noting that as already pointed out by 
Peierls for phonon-phonon interactions the major physical difference is between one dimension and higher. Three- 
phonon interactions in one dimension are severely restricted and usually very weak). However, the derivation can be 
easily generalized to three dimensional systems. Because it is not a major goal of this work to achieve quantitative 
agreement between theoretical results obtained with this phenomenological model and experimental data, we restrict 
the three-dimensional consideration to the same vibrational modes (one acoustic and one optic dispersion sheet). 
Note also that for finite temperature (T > To) three wave broadening, the integrals entering the self-consistent 
equations for the broadenings, only weakly depend on the region of wave vectors in the vicinity of the boundary of the 
Brillouin zone. Therefore, aiming only at a qualitative description of experimental data, in the first approximation 
one can perform this integration over an appropriately (i.e., with equal volumes) chosen spherical volume. Combining 
everything together we end up with the following three dimensional self-consistent equations for the mode broadening 
in the guest-host system. 



Ag2 



"'1(9) 



1 /.tt/2 



71(9) = — H / dt / dpp^/i^(p) 



1 Jo 
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/ /i(g)(7i(p) +7i(g))K(g) - wi(p)) 
\{wi{q) +wi{p) - wi(s))2 + (71 (p) +7i(s))2 

, /f(s)(7i(p)+7i(g))K(^) + ^i(p))/2 
(u)i((7) - - u;i(s))2 + (71 (p) +7i(s))2 

/l(g)(7i(p)+72(^))K(g)-i^i(p)) 
{wi{q) + wi{p) - W2{s))'^ + (71 (p) + 72(s))2 

72(g) - ^ / dt / dppVf (P) 



u'2(g) j-i JO 

■ {w2{q) + wi{p) -W2{s))'^ + (71 (p) +72(5))^ 

/f(-'^)(7i(p)+7i(^))(^'i(^) + ^i(p))/2 
{W2(q) - wi(p) - wi(s))2 + (71 (p) +7i(s))2 

/l(g)(7i(p) + 72(g) )(W2(S) +wi(p)) 
{W2iq) - wi{p) - W2(s))2 + (71 (p) +72(s))2 



(10) 



where s — \/ cp- + p^ — 2qpt. 



The self-consistency equations ([TOl) include as pointed out above only two modes: one optic and one acoustic (the 
latter one is longitudinal by its construction). Thus, we have in particular neglected two transverse acoustic modes. 
Lax et al., [23i] have shown that a given acoustic phonon cannot decay into other modes with higher velocity at any 
order in the anharmonicity. For the lowest-lying acoustic phonon mode one then expects anomalously long lifetimes. 
Similarly the phonon Umklapp processes have negligible effect on the heat flux carried by the low-energy phonons 
with small wave vectors. The latter would lead to very high thermal conductivity. For this reason, calculating the 
phonon broadening and thermal conductivity, we cut-off corresponding integrals (see below) at a certain small wave 
vector. Being interested in the only relative effects upon inserting guest atoms into the host lattice, this, and some 
other approximations made above driven by pure desire to make formulae simpler and to clarify physical ideas, do not 
affect our conclusions on the thermoconductivity reduction qualitatively. Moreover all these assumptions are easy to 
relax, if needed, for example other phonon-scattering mechanisms, providing finite broadening for long-wave phonons, 
can be considered on the same footing. 



IV. NUMERICAL RESULTS 



The self-consistent integral equations (|10p are our main results and they are ready for further inspection to find 
their solutions. Solving the system pH]) we find the broadening of the spectrum. Luckily for us the solutions to these 
rather complicated looking equations (e.g., they are not symmetric with respect to 1 ^ 2 exchange, as one could 
expect since there is no symmetry relation between acoustic and optic modes) can be found quite fast numerically (a 
few minutes for PC even using standard Mathematica software). Thus we will present results in the form of computed 
Mathematica plots. 

Figure O presents the results for the following parameter values kq = ki = 2, X = 0.8. Each branch is presented by 
three lines, which correspond to w{q) and w{q) ±7(5). Figure [S] shows the calculated (arbitrary normalized) density 
of states for the same parameter values. We should add that in the low-energy (frequency) region the vibrational 
density of states has to follow the well-known Debye law oc w^. The wiggles shown in the Fig. [6] on top of such a 
parabolic function are due to numerical artifacts. They can be easily eliminated by merely including more points in 
the integration code. Since we are not interested in precise absolute numbers but only in the more or less correct 
shape of the curves, we restrict ourselves to only rather crude numerical calculations. In the harmonic approximation, 
the peaks in the vibrational density of states (for the parameters corresponding to Fig. [6]at ~ 1.5, and llJ2 — 2.5) 
transform into Van Hove singularities (see the Fig. [7]). Figure [5] presents the spectra for the parameter values kq — 2, 
Ki = 1, and the same triple vertex A — 0.8. It is easy to see, that the broadening of the spectrum in this situation 
where the guest vibrations have softened with respect to the previous case is significantly larger. Figure O shows 
the calculated density of states for the same parameter values. Figure [10] shows (similarly to Fig. [7]) the harmonic 
counterpart of the density of states shown in Fig. [5] We should add that the value of the triple vertex used in our 
calculation was chosen such as to produce a broadening similar to that observed experimental in the case of soft guest 
vibrations, i.e. such that for large wave-vector values the relative width of the phonon line "f/wi 1/8. 
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FIG. 5: Spectrum and broadening for the force constants no = 2, ki = 2. 
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FIG. 6: Density of states, corresponding to Fig. (5] 
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FIG. 7: Harmonic density of states (dashed lines). Solid line shows the effect of anharmonic smearing on the singularities. 
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FIG. 8: Spectrum and broadening for force constants kq = 2, ki = 1. 




FIG. 9: Density of states, corresponding to Fig. [8l 



With all characteristics of the vibrational spectrum in hands, the thermal conductivity of the system rj can be 
calculated by the following expression derived from the Boltzmann kinetic equation [isj . [iJI 



C ^ 



dqq^ 



dwj {q) 
dq 



1 



(11) 



where C is the specific heat of the solid and the sum is taken over the two branches of the spectrum. Unfortunately, 
bluntly applied, this in the literature widely used formula for the thermoconductivity rj, is wrong. Indeed, in the 
derivation of this expression it has been assumed that the three-phonon mechanism of thermal conductivity is op- 
erational for all phonon branches, because the energy fluxes carried by the various groups of phonons are additive, 
and the same is true for their contributions to rj. If that mechanism does not work for one group of phonons, it 
would not provide finite thermoconductivity. However, the expression (jlip yields to a finite rj even if there were 
only normal (conserving quasi-momenta) processes, in contradiction to the notorious Peierls theorem (l6j . claiming 
that normal processes yield to 77 = oo. Luckily in our case (isotropic three dimensional guest-host system with two 
modes with dispersion laws ([4]), this issue is not relevant. The predominant contributions into the integral ()11|) do not 
come from very small wave vectors (excluded by our cut-off procedure), and for wave vectors larger than this cut-off. 
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FIG. 10: Notations are the same as in Fig. [7] 



just the Umklapp processes hmit the thermal flux of higher energy phonons. The physical picture of heat transfer 
looks as follows. Normal phonon - phonon collisions, in which the total quasi-momentum is conserved establish only 
local internal equilibrium in the phonons, and the Umklapp processes adapt this local equilibrium to the temperature 
gradient, and establish a finite value for the thermal conductivity rj. With said above in mind we are calculating 
numerically rj accordingly with Exp. (jlip . 

We find that the thermal conductivity in the case, when guest vibrations are soft, turns out significantly smaller, than 
in the pure host system 77(^1 — ko/2) ~ Q.\rj{Ki — kq). This result is independent of the absolute value of the force 
constants and equally of the anharmonic coupling constant. Thus, based on our calculations we anticipate that even 
a relatively modest softening of the guest atom vibrations should yield a noticeable reduction of thermoconductivity. 
We do believe that understanding of the underlying mechanisms is essential to predict the behavior. 



V. CONCLUSIONS: NEW OPPORTUNITIES FOR STUDYING ? 



Summarizing, in this paper we have investigated certain peculiar features of vibrational dynamics and thermal 
transport in guest-host systems. Since these systems are rather complicated and depend on many not easily accessible 
microscopic parameters, any theoretical study has to rely first on very crude simplifications. The model investigated 
theoretically in our paper includes only the minimal ingredients necessary to capture generic features exhibited by 
the compounds studied by experiment. We consider only two kinds of particles: strongly coupled host atoms, which 
form the structural cavities, and guest atoms confined in these cavities. In spite of the confinement, for relatively 
small amplitude vibrations of the guest atoms, their interaction with the host particles is weaker than between the 
host atoms. 

One of the main difficulties in comparing the results of our effective and simplified model with specific experimental 
measurements or elaborated microscopic numerical simulations is the availability of an accurate connection between 
experimental control parameters and theoretical model coefficients. The actual values of the parameters are determined 
by the microscopic (atomic level) interactions, and are in fact unknown. Without prior knowledge of the actual values 
of microscopic parameters we follow the logic of simplification. In our model the whole physics is lumped into the 
effective spring constants kq and ki and the three phonon interaction parameter A. We calculate the vibrational mode 
(phonon) dispersion laws, and self-consistently the line broadenings for the empty and filled (by the guest atoms) 
cavities. We show that the soft phonon sheet (which appears for guest atoms weakly bound in the cavities) yields large 
mode broadening via three wave anharmonic mode coupling. In turn the Umklapp part of this broadening reduces 
drastically the phonon thermoconductivity coefficient. We conclude that for a system with a soft mode vibrational 
spectrum the broadening is always considerably larger (thermoconductivity is always considerably smaller) than for 
a system without this soft mode (provided all other characteristic are the same). 
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Many of the points made above can be found in the hterature, and approaches similar to our paper have been made 
by other authors (see e.g. the review paper [l[ and references therein). However, some important differences to our 
work should be noted. We focus our study (and investigate in some detail) on the interplay between softness of the 
guest mode vibration and the host atom elasticity. Existing theoretical treatment of thermal conductivity is based on 
molecular dynamics simulations 3] which should utilize model dependent microscopic interaction potentials, and the 
Green-Kubo relation for the extraction of the thermal conductivity from the heat current correlation functions. Our 
calculations are much simpler and rely on a model (with only two relevant phenomenological parameters) that is based 
on general and qualitative assumptions. We feel that it accommodates a picture of the guest-host system vibrations 
and thermal transport without going into finer details. Of course we are left with many questions unanswered which 
must be pursued in further work. More elaborated models including many particles in the elementary cell, a realistic 
(corresponding to the crystalline symmetry) Brillouin zone etc., can be treated within the same approach and with 
the same conceptual ingredients. For partially filled cavities in the gust-host systems, special care should be taken 
about defects (both, point-like and extended imperfections, like pores and grain-boundaries). Luckily some general 
conclusions can be obtained without calculations. First note that while the specific heat and group velocity entering 
the expression for the thermoconductivity ()11|) are not sensitive to crystal imperfections, the mean free path is. For 
point defects which are small compared to the wavelengths of interest, the phonon scattering cross section cx Auj^, 
(where the coefficient A is proportional to Vb(AA//M)^ with Vq effective volume of the defect, and AM is the mass 
difference between the defect, e.g., guest atom, and the host atom mass M). Therefore point defects effectively scatter 
high frequency phonons (however, it should be noted that this expression strictly applies to a random collection of 
dilute point defects, whereas the guest molecules are neither dilute nor randomly arranged). Because mainly waves 
of a high group velocity make a substantial contribution to thermal conductivity, in our case mainly low frequencies 
play a role. It means that point defects are not very essential for our analysis. Large size imperfections scatter 
independently of frequency with a scattering cross section of the order of their geometrical size. These separation of 
frequency scales would alow in principle to distinguish this contribution to the reduction of thermoconductivity from 
guest-host soft mode mechanism of the thermoconductivity reduction, investigated in our paper. In such a generalized 
model the precise mode spectrum turns out to be much more complicated than the one found in our paper. However, 
the additional subtleties do not modify the general picture beyond the thermoconductivity reduction. Fact is that 
the parameter that tunes the anharmonic coupling in the high temperature limit T > Tjj only scales the final results. 
Therefore we anticipate that in the optimum conditions (see section |lT] above) , the relative reduction of the thermal 
conductivity will be robust with respect to the made approximations and almost A (see the Eqs. (1101) ') independent. 
In this sense our minimal model is just at the border between those that are too primitive to fit even qualitatively 
the data, and those that fit the data too well by using too many parameters. Our model may be an appropriate tool 
to generate predictions that can be in principle tested experimentally. 
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